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We consider the massive Dirac neutrino electric charge and magnetic moment within the context of 
the standard model supplied with SU(2)-singiet right-handed neutrino in arbitrary gauge. Using 
the dimensional-regularization scheme we start with the calculations of the one-loop contributions to 
the neutrino electromagnetic vertex function exactly accounting for the neutrino mass. We examine 
the decomposition of the massive neutrino electromagnetic vertex function. It is found by means 
of direct calculations that the massive neutrino vertex function contains only the four form factors. 
Then we derive the closed integral expressions for different contributions to the neutrino electric form 
factor, electric charge, and magnetic moment. For several one-loop contributions to the neutrino 
charge and magnetic moment, that were calculated previously with mistakes by the other authors, 
we find the correct results. We show that the electric charge for the massive neutrino is a gauge 
independent and vanishing parameter in the first two orders of the expansion over the neutrino 
mass parameter b = (m„/Mw) 2 . From the obtained closed two-integral expression for a massive 
neutrino electric form factor it is also possible to derive the neutrino charge radius. In the particular 
choice of the 't Hooft-Feynman gauge we also demonstrate that the neutrino charge is zero in all 
orders of expansion over b, i.e. for arbitrary mass of neutrino. For each of the diagrams contributing 
to the neutrino magnetic moment, we obtain the expressions accounting for the leading (zeroth) 
and next-to-leading (first) order in b, where the gauge dependence is shown explicitly. Each of the 
contributions is finite and the sum of all contributions turns out to be gauge-independent. Our 
calculations also enable us to obtain the neutrino magnetic moment in theoretical models that differ 
from each other by the values of particles' masses, including the case of a very heavy neutrino. The 
general expression for the massive neutrino magnetic form factor is presented. 

PACS numbers: 13.40.Gp, 13.40.Dk, 14.60.St, 14.60.Pq 



I. INTRODUCTION 



The recent experimental studies of the astrophysical and terrestrial neutrino fluxes provide the convincing evidences 
for the non- vanishing neutrino mass and neutrino mixing pj . These properties of neutrino are attributes of the physics 
beyond a scope of the standard model. An important information on the structure of the future model of interaction 
can be obtained in the investigation of radiative corrections to the properties of neutrino that in principle can be also 
verified in experiments. A critical test of a theoretical model is provided by the direct calculation of such characteristics 
of neutrino as its electric charge and magnetic moment. In this respect it is interesting to examine the gauge and 
neutrino mass dependence of these quantities. 

The electric charge and magnetic moment are the most important static electromagnetic properties of a particle. 
Their values are determined by corresponding form factors if the external photon is on a mass shell. In spite of the 
fact that the electromagnetic form factors are not measurable properties of a particle at nonzero momentum transfer, 
there are processes where the off-shell external photons are important. The example is the radiative corrections to 
the fcrmion-fcrmion scattering. The fcrmion electromagnetic vertex function and, in particular, its representation in 
terms of form factors in the case of off-shell external photon has been considered in Refs. 0, 01 within various gauge 
theories. 

To the best of our knowledge, however, there is no direct one-loop calculations of the neutrino electromagnetic 
vertex as well as of the neutrino charge and magnetic moment which are performed within the context of the standard 
model in the renormalizable gauge and which explicitly take into account the neutrino mass. It is worth to be 
noted here that the massive Majorana neutrino cannot have neither magnetic nor electric dipole moment. Due to the 
lepton flavour non-conservation that has been confirmed in the neutrino experiments, a neutrino could have flavour- 
off-diagonal transition magnetic moment, which is also allowable for Majorana neutrino. The differences between the 
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Dirac and Majorana electromagnetic properties are explained in detail in [j]. 

The neutrino vertex function in the limit of small neutrino mass was considered in Ref. [5j- There are several 
works where the neutrino charge is calculated within the standard model using the unitary, linear R^, and 't Hooft- 
Feynman gauges 0, H> El ■ The one- loop contributions to the neutrino magnetic moment in the standard model 
have been also considered previously 0, 0, ITTl IT^] . The vanishing of the massless neutrino charge emerges from 
the unbroken electromagnetic gauge invariance. The corresponding Ward identity has been derived in Refs. ItH] 
using the background field method. In the recent studies |l4j a computation of one-loop electroweak diagrams, which 
contribute to the neutrino charge and magnetic moment in the background field method and in the linear R^ gauge, 
is presented. However, all the previous calculations of the neutrino charge have been performed under assumption of 
a vanishing neutrino mass. With respect to the neutrino magnetic moment, an according treatment of this quantity 
was also performed to the leading order in the neutrino mass that is valid for the case of neutrino being much lighter 
than the corresponding charged lepton, m Ue <C rag. In addition, the results presented in Ref. fl4| for the gauge-fixing 
parameter dependence of several one-loop contributions to the neutrino charge and magnetic moment are incorrect. 

In this paper we consider the massive Dirac neutrino charge and magnetic form factors in the context of the stan- 
dard model supplied with SU(2)-singlct right-handed neutrino. Using the dimensional-regularization scheme, we start 
(Section [Hj with the calculations of the one-loop Feynman diagrams that contribute to the neutrino electromagnetic 
vertex function A At (q r ) in the general R^ gauge. It should be noted that, contrary to the previous studies, we ex- 
plicitly account for non- vanishing neutrino mass. In Section III Al we examine the structure of the massive neutrino 
electromagnetic vertex function. The decomposition of a fermion vertex function in terms of the four well-known 
electromagnetic form factors (presented, for instance, in Refs. 0, 0) nas been established using general principles 
such as the Lorentz and CP invariance and the hermicity. We analyze this decomposition and verify it by means of 
direct calculations in the case of massive neutrino within the standard model supplied with SU(2)-singlet right-handed 
neutrino. Such direct calculations were never undertaken previously. 

We present the general expressions for the contributions to neutrino electric form factor in Section II I II Then in 
Section ITlI Al we study the neutrino electric charge and analyze neutrino mass and gauge dependence of corresponding 
contributions arising from different Feynman diagrams. Although there is no doubt that the neutrino electric charge 
within the standard model is a gauge independent and vanishing quantity, however this fact has not been yet actually 
demonstrated in the case of a massive neutrino. Our calculations allow us to determine the neutrino mass and 
gauge parameter dependence of the one-loop contributions to the neutrino charge. We also obtain a corre ct g auge 
dependence for the contributions of several diagrams to the neutrino charge that have been calculated in Ref. with 
mistakes. Within the one-loop level we show that the neutrino electric charge is gauge independent and vanishing in 
the "zeroth" and first order of the expansion over the neutrino mass parameter b — (m u /Mv/) 2 ■ Moreover, for the 
particular choice of the 't Hooft-Feynman gauge we also demonstrate that the neutrino charge is zero for arbitrary 
neutrino mass, i.e. in all orders of the expansion over the parameter b. The obtained formulae can be used for studying 
the massive neutrino charge radius fSection |III C|) . 

In Section IIVI we consider the neutrino magnetic form factor using the one-loop contributions to the neutrino 
electromagnetic vertex derived in Section ^ For each of the contributions to the neutrino magnetic moment we 
derive the integral representations that exactly account for the gauge-fixing parameters as well as for the neutrino 
mass and corresponding charged lepton mass parameters (b and a — (mi/M^) 2 ). Then for each of the diagrams 
we perform an integration and obtain the explicitly gauge-dependent contributions to the neutrino magnetic moment 
accounting for the leading (zeroth) and next-to-leading (first) order in the expansion over neutrino mass parameter 
b. The sum of all the contributions turns out to be gauge parameter independent. However, our results for several 
contributions to the neutrino magnetic moment in the leading order in the neutrino mass disagree with those of Ref. [LH ] 
for the corresponding contributions. In particular, contrary to the results of Ref. [T^| . not all the contributions are 
gauge independent. Our calculations enable one to reproduce the correct value for the neutrino magnetic moment 
in any gauge including also the unitary gauge for which the results of Ref. jl4| are incorrect. In this Section we get 
final expressions for the massive neutrino magnetic moment in the various ranges of neutrino, charged lepton and W 
boson masses: m„ <C mg <C Mw, m£ <C m v <C Myy, and m,£ <C Mw "C m v . The last case amounts to a very heavy 
neutrino that is not excluded by the LEP data (see, e.g., Ref. [lf|). We also discuss the general formulae for the 
massive neutrino magnetic form factor at non-zero momentum transfer. 

The conclusions are made in Section We also include a list of the Feynman rules (Appendix [SJ and the typical 
Feynman integrals (Appendix used in our calculations. 



II. VERTEX FUNCTION OF MASSIVE NEUTRINO 



The matrix element of the electromagnetic current between neutrino states can be presented in the form 

(v(pi)\J™Hp))=u( P >)A,( q )u(p), (2.1) 
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where the most general expression for the electromagnetic vertex function A M (g) reads 



A m(?) = /q(? 2 )7p + fM{q 2 )i<J^q" - fE(q 2 )cr^q u j 5 + fA{q 2 )(q 2 Jp, - q/j,A)l5- 



(2.2) 



Here fq{q 2 ), fuiq 2 ), /b(9 2 ) an d /a(<? 2 ) are respectively the electric, dipole electric, dipole magnetic, and anapole 
neutrino form factors, q^ = p'^ — p^, = (i/2)[y fMl 75 = — ij ^ 1 ^ 2 ^ 3 . Their values at q 2 — determine the 
static electromagnetic properties of the neutrino. In the case of Dirac neutrinos, which is considered in this paper, the 
assumption of CP invariance combined with the hermicity of the electromagnetic current J^ M implies that the electric 
dipole form factor vanishes. At zero momentum transfer only /q(0) and /a/(0), which are called the electric charge 
and the magnetic moment, respectively, contribute to the Hamiltonian H m t ~ J^ M A^ that describes the neutrino 
interaction with external electromagnetic field A^. 

There is an important difference between the electromagnetic vertex function representations in the cases of massive 
and massless neutrino, respectively. If we consider a massless particle, from Eq. I|2.2|l it follows that the matrix element 
of electromagnetic current can be expressed in terms of only one form factor (see, for example, Ref. |l6j |) 



u{p')^n(q)u(p) = /i)(g 2 )M(p')7M( 1 + 75)u(p). 



(2.3) 



Thus, the electric charge and anapole form factors are related to the function /d(<7 2 ) by the trivial identities 



f Q {q 2 ) = f D (q 2 ), f A {q 2 ) = f D (q 2 )/q 2 



(2.4) 



However in the case of a massive particle, there is no such simple relation between the electric and anapole form 
factors since we cannot neglect the q^ ^75-matrix term in the anapole form factor. Moreover, the direct calculation 
of the massive neutrino electromagnetic form factors shows that, besides the ordinary electric charge and magnetic 
moment, each of the Feynman diagrams gives non-zero contribution to the term proportional to 7 M 7 5 -matrix. These 
contributions does not vanish even at q 2 = 0. This problem is related to the decomposition of the massive neutrino 
electromagnetic vertex function. Taking into account the importance of this problem, we present the direct calculation 
that verifies the decomposition given by Eq. (|2.2(l . Using the developed in the next section technique for studying 
the neutrino electric charge, we find that the sum of contributions of the complete set of Feynman diagrams to this 
additional additional "form factor" is zero at q 2 = 0. The vanishing of the considered "form factor" at q 2 ^ for the 
particular choice of the gauge is also demonstrated in the next subsection of this paper. 

We present below the one-loop calculation of the electric charge and magnetic moment of the massive neutrino within 
the context of the standard model supplied with SU(2)-singlet right-handed neutrino in the general gauge. The 
one-loop contributions to the neutrino electromagnetic vertex A M (q) are given by the two types of Feynman diagrams: 
the proper vertices (Fig. l(a)|l(f)' I and the 7 — Z self-energy diagrams (Fig. 2(a)|2(hj| . We use the Feynman rules 
given in Appendix^to find the contributions to the neutrino vertex function A^(q). In the dimensional- regularization 
scheme the contributions of the proper vertices diagrams (Fig. l(a)|l(f) I can be written as 
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FIG. 1: |(a)|(f)| proper vertices. 
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where m„, and are the masses of neutrino, W boson, and the charged isodoublet partner of the neutrino, 
respectively, e is the proton charge, g is the coupling constant of the standard model, 9w is the Weinberg angle, 
a = l/£ is the gauge parameter of W boson, Pl,r = (1 ± 7s)/2 are the projection operators. 



The contributions of 7 — Z self-energy diagrams (Fig. 2(a)|2(h) I to the vertex A M (g) are given by 
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FIG. 2: |(a)|(h)| The 7 - Z self-energy diagram. / denotes the electron, muon, and r-lepton as well as u, c ,t, d, s, and b quarks. 
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Here Mz and az are respectively the mass and gauge parameter of Z boson. In Eq. I|2.17f] . "— " and "+" stand for 
the "upper" (u, c, and t quarks) and "lower" (electron, muon, T-lepton as well as d, s, and b quarks) components of 
an isodoublet, nif and Q f are the mass and electric charge (in the units of e) of a fermion circulating within the loop. 

It is convenient to decompose each of the 7 — Z self-energy contributions at arbitrary q 2 and explicitly extract the 
transversal term: 
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Using Eqs. 12.11fl - 12.17fl for the contributions of the 7 — Z self-energy diagrams in the form of the Feynman integrals 
as well as Eq. (I2.18fl . it is possible to present the functions A^(a, q 2 ) and B^\a, q 2 ) (j = 7, . . . , 14) in the explicit 
form 
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In the derivation of Eqs. H2.19|) - l|2.32|) we have used the properties of the 7 matrix algebra in N dimensional space 
and the expressions for the characteristic loop integrals presented in Appendix FBI 



8 



A. The decomposition of massive neutrino electromagnetic vertex function 



In the direct calculations of the massive neutrino electromagnetic vertex function, taking into accounting the 
complete set of the Feynman diagrams, one reveals that, besides the well-known four terms, in Eq. (|2.2[1 appears an 
additional term proportional to 7p7 5 -matrix. Therefore we introduce the additional "form factor" /s((? 2 ). In this 
subsection we analyze this "form factor" and show by explicit calculations that f${q 2 ) = for arbitrary q 2 and for 
particular choices of the particles' masses and gauge-fixing parameters. 

Let us first consider the val ue of /s(g 2 ) at q 2 = 0: ip — f<i{q 2 = 0). The contributions to the "charge" ip of the 
proper vertices diagrams (Fig. l(a)|l(f) I have the form: 
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are the charged lepton and neutrino mass parameters, respectively, D a = a + (a — a)z — bz(X — z) and D = D a =\ = 
a + (1 — a)z — bz{\ — z). 

In Eqs. (|2.33() - (|2.37() we assume the mass parameters a and b as well as the gauge parameter a to be arbitrary. 
We have calculated the integrals in Eqs. H2.33fl - H2.37fl . however the results, being expressed in elementary functions, 
are rather cumbersome. Therefore we also perform corresponding integrations in the first two terms of expansion 
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over the neutrino mass parameter b for arbitrary values of the charged lepton mass parameter a and the gauge-fixing 
parameter a. In this case the sum of the proper vertices diagrams contributed to the "charge" ip can be written as 



^prop.vert.)^^) = £^=M^ (a, a) + btpf (a, a) + <D(b 2 )}. (:>.:Wi 



i=l 



The contributions of the 7 — Z self-energy diagrams (F ig. |2(a)|2(h) 1 to the " charge" tp coincides with those to the 
neutrino electric charge Q^~ z ^> and thus are given by Eq. I|3.26|l ( the details of the neutrino electric charge calculations 



can be found in Section lill A|) . We have calculated the functions <p£ (a, a) and <pf (a, a) and found out that the sum 

of all the contributions (p^ (a, a) exactly cancels the contribution of the 7 — Z self-energy diagrams. This our result 
corresponds to the case of a massless particle. Therefore the "charge" <p of a massless neutrino is zero. Then, summing 
the contributions (pf' '(a, a) we reveal that the value of the "charge" <p is also zero in the next order of the expansion 
over the neutrino mass parameter b. 

Now let us consider the value of the "form factor" f${q 2 ) at non-zero momentum transfer. In the subsequent 
calculations we have to fix the gauge in order to simplify the formulae. We set otz — 00 and a = 1 that corresponds 
to the unitary gauge for Z boson and the 't Hooft-Feynman gauge for W boson. However, we do not restrict 
ourselves considering either light neutrino or light charged lepton: the mass parameters a and b are arbitrary in all 
our calculations. In this case the function B(q 2 ) in the decomposition of the 7 — Z self-energy diagrams [see Eq. <|2.18ll ] 
takes on the form 

14 

B{q 2 ) = J2 bU) (v 2 ) = 2cos<V sin6W M^M 2 z G f {-2uj + g c (r) cos 2 9 W + g s {r) sin 2 6 W ), (2.40) 

3 = 7 



where 



7 It 



[ dx{6x 2 -&x + l) ln(l - C) - 9 [ dx ln(l - C), 
Jo Jo 



T T I 1 . , 



5 (t) = - + -/ dx{6x 2 -6x + l)ln{l-0- I dxln(l-C). 



(i 



We note that g s (0) = g c (0) = 0. To show that <7 s (t) = 3c(t) for any value of r we consider their difference 
g(r) — .9s(t) — g c {T~)- The function g(r) can be represented as follows 

g(T) = ^T + 8[ cfom(l-C)+4T / da;(6a; 2 -6a; + l)ln(l-C). 
3 Jo Jo 

Expanding ln(l — Qj in the formal series 



00 u 



k=0 



and performing the integration with the help of the formula 



{l + s + 2) 



we receive that g(r) = for any value of r. Thus, the expression for B(q 2 ) can be represented in the form 

B(q 2 ) = 2 cos 6 W sin ewM^M^Gp (-2w + g s (r)). 



It is interesting to note that, in contrast to Eq. (|2.4UI) . in this expression the dependence on cos 9w an d sin 9 W 
is absent. We have also verified that the analogous property of the function B(q 2 ) remains valid in the case of an 
arbitrary gauge. 

The 7 — Z self-energy diagrams contribute to the fsiq 2 ) "form factor" according to the formula 

h (q) ~ 4Af 2 cos^ M ^ <Z >■ 



10 



Now we turn to the contributions of the proper vertices diagrams (shown in Fig. l(a)|l(f)' I to the /^(q 2 ) "form 
factor" at arbitrary q 2 : 



(prop. vert. 



(9 2 ) = £/ 5 ( V) 



4tt 2 V2' 
a — b 



2 

a — b 



M^i-uj+f dz f dy\n(D' - ry(z - y)) - 3 f dz f dy\n{D-Ty{z-y))+ 
I Jo Jo Jo Jo 

dz dy\n{D' -Ty{z~y))~ dz dy\n{D-Ty(z-y)) + ^j + 



dz / dy— 



1 



D' - ry(z - y) 
1 f z 

dz / dy 



-(a- 6(1 - z)') + r ( -z{\ - z) + y(z - y) 



1 



D - ry(z - y) 



1 



1 



1 



a-bz +t \3y(z-y)- -z- -z + ~y - ~(a - b)z +{a-b)y{z-y) 



(2.42) 



where D' — I + (a — l)z — bz(l — z). Eq. (|2.42|) can be analyzed in the same manner as we have treated the function 
g(r). For instance, let us present the calculations of one of the integrals in Eq. (|2.42() 



dz / dy(a — bz 2 ) 



I(r)- 

We again expand the integrand in Eq. I|2.43|l in the formal series 

1 



D - ry(z - y) ' 



(2.43) 



D - ry(z - y) 



1 I 



fc=0 



Then, we carry out the integration over the variable y using Eq. (|2.41l) . The obtained result should be transformed 
according to the identity 



1 z l 
dz Jv^i( a - bzZ ) = T + 



, 2 s 1 k — I — 1 



dz—-, k > 1, I > 0, 

JJk ' — ' — ' 



D fc + lv "~ ' k ' k 

which can be proven by means of partial integration. Finally, we get the following expression for the function I(r): 

- k (fc!) 2 



IM = V — — W - V fc + 2 ^ r k 
y 1 ^ fc (2fc + l)! fc (2fc + l)! J 



1 z 2k+l 

dz- 



D 



(2.44) 



Note that the first term in Eq. 12.44JI does not depend on the charged lepton and neutrino mass parameters, a and b. 
It is this term that cancels the corresponding contribution of the 7 — Z self-energy diagrams. The subsequent analysis 
of the remaining contributions of the proper vertices diagrams can be performed in the same manner as we have done 
it for the function I(t). Finally, we obtain that 

/ 5 (g 2 ) = / 5 (7 - z) (9 2 ) + / 5 (prop - vert - ) (? 2 ) = o 

for any value of q 2 and for arbitrary charged lepton a and neutrino b mass parameters. It should be noted that 
the decomposition of the fermion electromagnetic vertex function in terms of the four form factors was established 
previously with the use of only general principles such as the Lorentz invariance and hermicity of the electromagnetic 
current operator. We have demonstrated the validity of this decomposition by means of the direct calculations of the 
corresponding Feynman diagrams. 



III. NEUTRINO ELECTRIC FORM FACTOR 



In this section we study the massive neutrino electric form factor. Using the results of the previous section for 
different contributions to the neutrino vertex A^(q) we extract in Eqs. I|2.5|l - (|2.17|l the coefficients proportional to 
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7 M -matrix that are, according to the decomposition l|2.2[l . the corresponding contributions to the neutrino electric 

form factor fq{q 2 )- 

First of all we consider the contributions of the one-loop proper vertices [Fig. l(a)|l(f)] to the neutrino electric 
form factor. Using the well known mass shell identity 

w(p')(Pm +Pn) u (p) = u(p')(2m„j fl - ia^ u q v )u{p), 

and carrying out an integration over the virtual momenta within the dimensional-regularization scheme (see for more 
details Appendix^ we derive the exact expressions for the contributions from the considered diagrams to the massive 
neutrino electric form factor in terms of the definite integrals 



^prop. vert.) _ 



z— 1 



where 



jti) , 9i cx, , 1 — a 



[ dz [ dymlDi- [ dz [ dy(a + b(l-z) 2 +T(l-z + y(z-y))) 7 ^-+ 
Jo Jo Jo Jo ®1 



dz / dy(bz 2 (a + 6(1 — z) 2 ) + ary(z — y)+ 
2 Jo Jo 

br(2zy(z - y)(l - z) + by{z - y) - z 2 (l - z)) + T 2 y(z - y)(l - z + yz- y 2 )) 



1 



1 



2>i(a) 23 1 



dz I dy(a + b + 6bz(l-z)+T(l-3z + 6y(z-y)))[lnS) 1 (a)-\n® 1 ], (3.1) 
2 Jo Jo 



/SV) 



a + b ( UJ 1 
2 + 2 



f dz / dylnSi(a) 
o Jo 



1 f 1 f z 1 

dz / dy(a 2 + abz 2 + b 2 z 2 - Aabz + ab + (a + b)ry(z - y))— -r -? , (3.2) 
o Jo 2Di(a) 



/qV) = \-\ - I dz I dylnS 2 (a)) +b j dz j dy(3az-az 2 -2a + bz(l-z)) 



o Jo 



o Jo 



S 2 (a) 



(3.3) 



1 

2)7 



/qV) = -w 7 (l + a)-l-3 / / dyln£ 2 +/ dz / dy(3fe(l - z) - r(z - y(* - y))) 
4 J J J J 

I f dz f dj/[(S 2 (a) + y(l-a))ln(S 2 (a) + y(l-a)) -S 2 lnS 2 ]- 
* Jo Jo 

dz / dy(2fo 2 (l - z) 2 (z(l - z) - y) - br(y(z - y)(5z - 3z 2 - 3y) + z(l - z) 2 - y(2 - y - y 2 ))- 
Jo 

r2 2/(^ ~ 2/)(l ~ z + 2/ z + 2/ + y 2 )) 



iD 2 (a)+tf(l-a)) Z»2 



for 



dz / dj/(36(l - z 2 ) + r(4 - 6(z - y) + lly(z - y))) [ln(S 2 (a) + y(l - a)) - lnD 2 ] 
1 Jo Jo 

1 1 2 



, dz / dy(bz{\ - 3z + z 2 + z 3 ) - ry(z - y)(z + z 2 - 2y)) 
2 Jo Jo 



£> 2 ©2(a) £> 2 (a) +y(l - a) 
dz / dy(6(9-13z + 4z 2 )-2ry(z-y)) [lnS 2 + lnS 2 (a) - 21n(23 2 (a) + y(l - a))] + 



T 

4./,, ./„ 

3r f 1 



dz / dy[2) 2 ln£) 2 +£l 2 ( a )ln2) 2 (a)-2(D 2 (a)+y(l-a))ln(D 2 (Q ! )+y(l-Q ; ))], (3.4) 
4 Jo Jo 
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(5) + (6) 



(r 



dy(a — bz) 



S) 2 (o)+»(l-a) 



rfz 



Here 



2)i(a) 
©2(a) 



a + (a — a)z — bz(l — z) + ry[z — y), 
a + (a — a)z — bz(l — z) + Ty(z — y), 



S5 2 (a)+2/(l-a) £ 2 ( a ). 
eZy(a + 5fc - 66z) [ln(D 2 (a) + y(l - a)) - ln£> 2 (a)] ■ (3.5) 



Si = 2)i(a = 1) = 1 + (a - l)z - - z) + ry(z - y), 
D 2 = £> 2 (a = 1) = a + (1 - a)z - 6z(l - z) + ry(z - y). 



Note that the values of the mass parameters of the charged lepton (a) and neutrino (b) are taken into account explicitly 
in Eqs. I|3.1[l - i|3.5[l . The gauge parameter a and q 2 are arbitrary in these formulae. 

The contributions of the 7 — Z self-energy diagrams (Fig. 2(a)|2(h) I to the electric form factor can be obtained 
using Eqs. I|2.10|l and (|2.18|) . Thus, one obtains 



g A^(a,q 2 ) + B^(a,q 2 ) 



4cos( 



'ir 



q 2 -M 2 



3 = 7, 



14. 



(3.6) 



Using explicit form of the functions A {l \a,q 2 ) [Eqs. l(2~TU|l - (|2~23|l ] and B^(a,q 2 ) [Eqs. g^St - jg^ ] as well as 
Eq. (|3.6I) . one can also derive the expressions for the contributions of the 7 — Z self-energy diagrams at arbitrary 
values of the gauge parameter a and q 2 7^ 0. 



A. Neutrino Electric Charge in Arbitrary Gauge 



In this section we consider the neutrino electric charge. At zero momentum transfer the sum of the contributions to 
the electric form factor determines the neutrino charge, /q(0) = Q. Our goal is to find its total value for the massive 
neutrino 

6 14 

z=l j=7 

and to study the mass (a and b) and gauge-f ixing (a a nd az) param eters dependence of the contributions from the 
different Feynman diagrams depicted in Figs. l(a)|l(f) 



and 



2(a)g(h) 

ma 



First we consider the one-loop contributions to the neutrino charge which arise from the proper vertices diagrams 
in Fig. l(a)|l(f)' Using of the more general formulae for the electric neutrino form factor Eqs. H3.lfl - 13.5fl . we obtain 
the exact expressions for the contributions from the considered diagrams to the massive neutrino charge in terms of 
the definite integrals: 



)W(a,6,a) = 



6Gf -_M 2 



4ir 2 V2 



a 

nl + 1 



1 — a 
12 



/ dz(l-z)hiD- 
Jo 



1 



dz(l - z)(a + bz 2 )— + - dz(l - zf{a + bz 2 ) 
i< D 2 J 



1 1 



dz(l — z)(a + b + 



6bz{l- z))[lnD a -\nD}+3 dz(l - z) [D a ha.D a - DlnD] j, (3.7) 



4tt 2 V2 



w 



( a + b foj 1 
l^(2 + 2 



dz(l — z) In D^j 

1 r 1 
2 



dz(l - z)(a z + abz 2 + b 2 z 2 - 4abz + ab) 



D a 



(3.8) 







j dz zlnDc^J 






Jo 



z (3az~az 2 -2a + bz(l- z))-^-\, (3.9) 
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Q (4) (a,6, a) 



eG F 



M^j-a;^(l + a)-l-3jf dzz\nD + 3b 



dz z 2 (l — z )~q~ 



dz / dy [{Da + y(l - a)) \n(D a + y(l - a)) - D In D] 



2b z I dz dy(l - z) 2 (z(l - z) - y) 



D a + y{l-a)) D\ 
ho J dz J dy(l-z 2 )[ln(^ Q +y(l-a))-ln^]|, (3.10) 



Q^+^(a,b,a) = ^%M^ 



4tt 2 V2 



dz dy(a~bz)— — ■ — 

o Jo D a + y(l-a) 



b I dz dy(l — z) (a — bz) 
lo Jo 



1 



D a +y(l-a) 

i J dz dy(a + 5b - 6bz) [ln(D a + y(l - a)) -In AJ j. (3.11) 



The integral expressions of Eqs. (|3.7|) - l|3.11|) for different proper vertices contributions to the neutrino charge exactly 
account for the charged lepton and neutrino mass parameters, a and b, and also for the gauge-fixing parameter a. 
We have calculated the integrals in Eqs. (|3.7|l - l|3.11|) . however the results, being expressed in elementary functions, 
are rather cumbersome. Therefore, we also perform corresponding integrations in the first two terms of expansion 
over the neutrino mass parameter b for arbitrary values of the charged lepton mass parameter a and the gauge-fixing 
parameter a. In this case the sum of the proper vertices diagrams to the neutrino electric charge can be written as 



Q (prop.vert.) (a)6)a) = ^=M^ (a, a) + bq® (a, a) + 0(b 2 )} 



(3.12) 



For 5q (a, a) we obtain 



1 



■ ( - 3a 2 + 4a 2 a In a - 5a 2 a + 2a 3 + 3a 3 - 3a 3 In a - 2a 3 In a- 



2 4(1 -a) 2 (a- a) 2 

2 2332122 23 3 3 4 

a a In a + 2a a — 3a a + aa + 6a a + 6aa — 3a — 4a a — 2a a + 4a a In a — 2a alna— 

2a 3 a 2 In a + 4aa 3 In a — 6a 2 a 2 In a + a 2 a 3 In a — 2a 2 a 2 In a + 3a 3 a 2 In a), (3.13) 



On 2 "* (a, a) = lo — I - — (2a 2 In a + 4aa — 4aa In a — a 2 + 2a 2 In a — 3a 2 

u 4 8(a - a) 2 v 



(3.14) 



(3) / \ a 



(a — a) ; 



■ (2a 2 In a + 4aa — 4aa In a — a 2 + 2a 2 In a — 3a 2 ) 



(3.15) 



% K a ) 



3 1 

""V +a > 8(l-a) 2 (a-a)(l-a) 



a 2 a — 4a 2 — 6a 3 In a — 6a 3 In a- 



lla 2 a + 5a 3 + 5a 3 + 5a 2 a 3 - 5a 3 a 2 + 10a 2 a 2 + lOaa - 6a 2 - 10a 3 a- 



6a a In a — 12a a In a + 12a alna + 12a alna + 6a a In a + a — a), (3.16) 



^ 5)+( Va) 



- (2aa In a + a 2 a — a 3 + 4aa In a — 



a 2 a 2 — aa + a 2 + a 3 a— 



4(a- a) 2 (l - a)(l - a) 

3a 2 In a + a 2 a 2 In a + 2a 2 a In a + 4a 2 a In a — a 2 a 2 In a — 4aa 2 In a — 2a 2 a In a — 3a 2 In a). (3.17) 
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Each of the coefficients q$(a, a), if considered separately, depends on the gauge-fixing parameter a and all of them 
[except for (a, a) and qjp(a, a)] are divergent. Note that, according to the expansion given by Eq. (|3.12|) . the sum 

6 



determines the proper vertex contribution to the charge in the limit of the massless neutrino. 

It should also be noted that the two diagrams of Fig. |l(e)|a nd 1(f) are convergent for every value of the gauge-fixing 
parameter a. One can check this statement using Ea. (|2.9)l . Indeed, these diagrams have the superficial degree of 
divergence equal to —1 and hence converge |17| . Therefore the according contributions to the electric charge must 
be finite as it is also shown by Eqs. (|3.11|) and (|3.17|) . Here we find a discrepancy with the corresponding results of 
Ref. 01 where the massless neutrino charge is calculated and the contributions of these two diagrams contain the 
ultraviolet divergencies. 

The next order over the neutrino mass parameter b of the proper vertex diagrams' contributions to the neutrino 
charge can be obtained if one expands the integrands in Eqs. I|3.7|) - (I3.11|I . keeps the terms proportional to &, and then 
carries out the integration. Taking into account that the functions D and D a also depend on 6, we find that 



£gf (o,a) = 0. 



(3.18) 



Thus, due to the fact that, as it is shown below in this section, the 7 — Z self-energy contribution does not depend 
on the neutrino mass, it follows that the neutrino charge term proportional to the neutrino mass parameter b is zero. 
Now let us turn to the 7 — Z self-energy contributions to the neutrino electric charge. The corresponding Feynman 



diagrams are depicted in Fig. 2(a)|2(h) Using Eq. I|2.18[) . which presents the decomposition of the functions Tlfll (q), 
as well as the explicit form of the functions A^\a,q 2 ) [Eqs. (|2.19|l - (|2.25|l ]. we find that 



0) = 0. 



Therefore only the terms proportional to B^\a, 0) are responsible for the neutrino electric charge in the 7 — Z 
self-energy Feynman diagrams and we have 



14 



Q 



(a) Q W) («) = - 



<) 



14 



3=7 



AMI cos 6 W 



3=7 



For each of the contributions Q^'(a) (j = 7 
n (7), s eG F 



14) from Eq. l|3I>|) we obtain 



4tt 2 v / 2 



Mw cos 2 6w{ u 



3+ -a(l + a) 
4 



5a 
~ ~8 



5a 2 



3a 3 In a 
~ (1 - a) 



(3.19) 



)(«) 



(a) 



3 + a 5 + a a /„ , a\ met 
~ 8^ 2 



( 1+ 2) (T^))' 



(3.20) 



g(9)( a ) = eGp Mir (cos 2 9 W - sin 2 Q w )-\ {-too. + a - a\ri a} , 
47r 2 V2 



(3.21) 



eGp 



(a) = -j-Jj= M 2 W cos 2 9 W \ - 7W (3 + a 2 ) + - + 



3 5a 2 3 



— -or In a 



(3.22) 



)(11) + (12)/„a _ eG F „, r2 „„„2 - 1 



'(«) 



~M w cos Oy/n {— + a — a In a} , 
47r 2 v2 2 



(3.23) 



Q {13) (a) = eGr M^(sm 2 6 W - cos 2 %)-{-wa + a-alna}, 
47r 2 V2 ' 2 



(3.24) 
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Q (14) («) = 0. (3.25) 

It is worth to mention that each of the contributions Q^(a) turns out to be independent on the neutrino (m v ) and 
charged lepton (me) masses. There is also no dependence on the masses m/ of the virtual fermions that circulate 
in the 7 — Z self-energy diagrams because of the properties of the 7-matrix algebra specified in Appendix [5] The 
dependence on the gauge- fixing parameter otz also cancels out within each of the contributions. Note that prior 
the integrations in Eqs. i|2.11fl - H2.17fl 1S carried out the az dependence drops out of each of the electric form factor 
contributions to the vertex function at arbitrary momentum transfer q 2 . 

Finally, for the sum of all 7 — Z self-energy contributions to the neutrino electric charge we have 

^<-, 7\, s &Gf . ,? f3 + a 5 + a alna / a\l . 
Q(l ^ = iWl M - {—» - — - 2(1=70 (' + 2 ) } ■ (3 ' 26) 

Some remarks should be made with respect to the divergent parts in Eqs. H3.13JI ~ H3.16JI anc - <|3.26J1 . The sum of all the 
coefficients in uj terms is zero, i.e. the electric charge of a massive neutrino vanishes for every number of dimensions 
N. The same property of the electric charge of a massless neutrino was determined in Ref. [9j. 

Now we can complete investigation of the neutrino charge in the zeroth order of the expansion over the neutrino 
mass parameter b summing together the contributions from the proper vertices given by Eqs. (|3.12JI - (|3.17Ji . 

Q(pro P .v CT t.) Mja) = ^L^ M 2 r J2q^(a,a), (3.27) 

v i— 1 



and Q^~ z \a) given by Eq. 

As a result we obtain that the neutrino electric charge in the zeroth order in the neutrino mass vanishes for every 
gauge in agreement with the final results of Refs. 0,0], where the calculations of the neutrino electric charge have 
been performed in the limit of vanishing neutrino mass. 

B. Neutrino Electric Charge in the 't Hooft-Feynman Gauge 

Within the 't Hooft-Feynman gauge it is possible to show explicitly that at the one-loop level the neutrino electric 
charge is zero for arbitrary mass of neutrino. The gauge-fixing parameter a = 1 in this gauge. Summing up the 
contributions of all relevant diagrams [Eqs. (|3.7Ji - H3.11JI and (|3.26J) ] we receive the exact expression for the neutrino 
charge at arbitrary values of the charged lepton and neutrino mass parameters, a and b, while a = 1: 

Q(a,b,a=l) = £^ M w{lJ Q dz[a(-(a + 2) + z(a + 4))+ 

ab(-l + z + z 2 - z 3 ) + 2bz 2 (l - 2z) + b 2 z 2 (l - z)]^+ 

J dz(l-Az)\nD+^-^ {\ + J ~ 2z ) ln -°) }• ( 3 ' 28 ) 

To analyze Eq. we use the formulae 

/ dz\a.D = -l+ ( dz(a-bz 2 ) — , [ dz zlnD = -- + - [ dzz(a~bz 2 ) — , (3.29) 

that can be proven by means of partial integration. Substituting Eq. 13.29JI into Eq. (|3.28JI we obtain that the neutrino 
electric charge vanishes for arbitrary neutrino mass in the considered gauge. 

C. Neutrino Charge Radius 

Using the closed expressions for the contributions to the neutrino electric form factor obtained above in this Section, 
it is possible also to derive the neutrino charge radius. Accounting for the next-to-leading term in the g 2 -expansion 
of the contributions in Eqs. (J3.1|) - IJ3.5JI and (J3.6|) . 



" 3 2 ||(0} 
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one can obtain the value of the massive neutrino charge radius as 

(O = -6©(0) 



and study its dependencies on the gauge and mass parameters. Here we should like to note that the problem of the 
massless neutrino charge radius has been discussed in details in 0, 0] . 



IV. NEUTRINO MAGNETIC MOMENT 



According to the general decomposition of the neutrino electromagnetic vertex function A M (g) given by Eq. (12.211 , 
the neutrino dipole magnetic form factor /m(<? 2 ) is the coefficient in the term proportional to ia liV q v . In this section 
we first determine /m(<Z 2 ) and then calculate at q 2 — the neutrino magnetic moment accounting for the two mass 
parameters (a and b) and for the gauge-fixing parameter (a) as well, 



fj,(a,b,a) = f M {q 2 = 0). 



(4.1) 



Note that the Feynman diagrams in Fig. 2(a)|2(h) do not contribute to the neutrino magnetic moment. Thus, the 
total one-loop value for the neutrino magnetic moment is given by 



/i(a, 6, a) = /j,^ (a, 6, a) 

i=l 



(4.2) 



la 



where fv-*' (a, 6, a) are the contributions to the magnetic moment from the corresponding diagrams shown in Fig 

EH 

We treat the neutrino magnetic moment in the similar way as we have analyzed the neutrino electric charge. Using 
Eqs. I|2.5|) - (|2.9|l for each of the contributions to the neutrino magnetic moment we receive 



jir 1 ^ (a, 6, a) 



eG F / f 1 , 2 1 1 ' ' 

=ra„< / dzz\\ — z 

4^ 2 V2 \ Jo D 2 



dz(l - z) 3 (a - bz) 



1 1 

15 a ~ 15 



ijT dz(l-z)(l-3z)[lnL> a -lnD] j, (4.3) 



// 2 ^(a, 6, a) 



eG F 1 

zUl, 



-=,,v v — , dz(l — z){— 3az + az 2 + 2a — bz(l — z)) — 

47r 2 v2 2 J D c 



(4.4) 



pt (3) (a, b, a) 



eGp 

=9 



If 1 1 

i v — / dz z(~3az + az 2 + 2a — bz(l — z))—— , 
2 Jo D a 



(4.5) 



/i^ 4 -* (a, b, a) = 



eG F 

=? 

4tt 2 V2 



dz dy(l~zY(z(l-z)-2y) 
1 Jo Jo 



D a + 2/(1 - a) 



1 

D 



dz / dy(-2 + 9z- Az 2 - 6y) [ln(D a + y(l - a)) - InD] \, (4.6) 



^+^(a,b,a) 



eG F / f 1 f* 1 

-=ra v <, dz dyy— ; r- 

4^ 2 V2 \Jo Jo yy D a +y{l-a) 



+ 



1 



dz / dy(l — z) 2 (a — bz) 



1 



D a + y(l-a) 
\J o dz J o dy(2-3z)[ln(D a + y(l-a))-lnD a }y 



(4.7) 
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It should be noted that these formulae exactly account for dependencies on the neutrino and charged lepton mass 
parameters (a and b) and the gauge- fixing parameter (a). 

To proceed further with the analytical calculations we expand the contributions to the neutrino magnetic moment 
[Eqs. H4.3fl ~ H4.7fl ] over the neutrino mass parameter b and consider the first two terms. Then from Eq. (|4.2fi we obtain 

M (a, b, a) = ^7=rn v £{/2«(a, a) + bp,? (a, a) + 0(b 2 )}. (4.8) 

* 2 — 1 

For each of the coefficients p!^\a,a) we have found the exact expressions in terms of algebraic functions, however 

they are again rather cumbersome. Therefore, let us consider more compact expressions for Ji^ (a, a) that can be 
obtained in expansion over the charged lepton mass parameter a. Thus, accounting for the terms up to the second 

order in a we derive for the coefficients p!^{a,a) 

m, . 2 10 — 3a + 61na — 61na „. , , „. 

/4>,«) = 3 + ^ a + 0(a 2 ), (4.9) 

-(2)/ \ 5 + 3 In a - 3 In a „. 2 , , A m\ 
Hl\a,a) = — a + 0(a 2 ), (4.10) 



^ 3)(a ' a) = T^ + 0(a " ) ' (4 ' n) 



5 a 



(a). . 2 — 7a — 3a In a + 5a 2 9 — 12a + In a + 5a In a + 3a 2 ,. , , „. 

4 } M = ^— w — a + G{a 2 ), (4.12) 



, , 1 — a + alna 5 — 16a — a In a + 11a 2 — 5a 2 In a , , v , 

4 5)+(6 W) = ^TT^P a + °^ ^ 



Eqs. If4.9fl - (f4.13fl together with Eq. If 4. 8(1 yield a value of the magnetic moment in the limit b — * that corresponds 
to the case of a light neutrino. We may compare our p,Q'(a,a) calculations with the results of Ref. ^4|. Our results 

for the contributions ^^(a, a), ^\a, a), and ^\a, a) disagree with those of the above cited paper. The Feynman 
diagrams corresponding to the contributions i = 5, 6 contain the unphysical charged scalar boson. This boson 
contributions should disappear in the unitary gauge when the gauge parameter a = oo. Thus, the contributions 
to the magnetic moment from these two diagrams must vanish in the limit a — > oo. This is exactly what we get 
from Eq. H4.13fl . However, the similar expression from Ref. [14| does not depend on the gauge parameter at all. An 
argument in favor of our results can be also obtained if one considers the value of the neutrino magnetic moment 
within the unitary gauge. Indeed, it is easy to show that using the results of Ref. 0] it is not possible to get the 
right value for the neutrino magnetic moment within this gauge. In the unitary gauge, only the diagrams shown in 



Fig. 1(a) and 1(d) contribute to the neutrino magnetic moment. The results for these two contributions, that can be 
obtained using the corresponding formulae presented in |l l| . are 



"» 4, ^-{b +o «° 2 »}' ( « 5) 

The sum of the leading terms in Eqs. If4.14fi and l|4.15fl differs from the well known result for the neutrino magnetic 
moment calculation (see, for example, Ref. 

3eGi?mj/ 
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This fact points out that the contributions of the three diagrams shown in Fig. l(d)|l(f) are calculated in with 



incorrect gauge-fixing parameter a dependence. Our calculation shows that to the leading order in m„ each of these 
three contributions of the diagrams in Fig. l(d)|l(fj are gauge- fixing parameter dependent. Note that calculations 



performed in |14j provide the correct results only within the 't Hooft-Fcynman gauge. 

Let us now consider the value of the neutrino magnetic moment in the "zeroth" order in the expansion over the 
neutrino mass parameter b taking into account all the contributions. The sum of the coefficients (14.9(1 - ()4. 13(1 is found 
to be independent on the gauge parameter a. The straightforward calculation of the neutrino magnetic moment in 
the limit 6 — > yields 

li (a,a) = F / - m v — l -rr(2-7a + 6a 2 -2a 2 lna-a 3 ) (4.16) 

47r 2 V2 4(1 - ay 

that is in agreement with Ref. [T^j . 

Considering the next-order over the neutrino mass parameter b contribution to the magnetic moment, we find out 
that the sum of the corresponding contributions of Eqs. ((4.3(1 - 1(4.7(1 to the coefficient /2i(a, a) is given by 

pi (a, a) = ^2fti (a, a) = — — — (5 - 26a + 6a In a - 36a 2 - 60a 2 In a + 58a 3 - 18a 3 In a -a 4 ). (4.17) 
i=i ^ a ' 

Thus, we explicitly show by Eqs. ((4.161) and l(4.17|l that in the one-loop level and to the second order in the expansion 
over the neutrino mass parameter b the neutrino magnetic moment is a gauge-independent quantity. 

The obtained Eqs. ((4.3|) - ((4.7|) also enable us to consider the magnetic moment of a rather heavy neutrino since the 
mass parameters a and b are arbitrary in these equations. Let the neutrino mass m v be much greater than the charged 
lepton mass mi (this case amounts to b 3> a). Approaching the limit a — > in Eqs. 1(4. 16(1 - 1(4. 17(1 . while keeping b 
constant, for the neutrino magnetic moment we receive 

3eG F f 5 , 1 , 

-.m v { 1+ — & + ••• \. (4.18) 



The recent LEP data require that the number of light neutrinos coupled to Z boson is exactly three. Any additional 
neutrinos must be heavier than 80GeV (see, e.g., Ref. [HI]). Using our formulas Eqs. ((4. 3|) - 1(4. 7(1 f° r a massive Dirac 
neutrino magnetic moment we can also examine the case of a very heavy neutrino. Let us consider the case of 
neutrino mass being even greater than W boson mass. To examine this situation we should fix the gauge parameter 
a in Eqs. 14.3|) - ((4.7I) for simplicity of the computations. In what follows we set a = 1 that corresponds to the 
't Hooft-Fcynman gauge. Thus, for the sum of all the contributions to the magnetic moment we obtain the expression 



2S w + 5t)^J dzz 2 ^, (4.19) 



V = z 2 - z(l - S w + S() + St, 

where we redefined our mass parameters and introduced the two new quantities: Sw = 1/b = {My/ /m y ) 2 and 
5 1 = a/b — {mi/rriv) 2 . The case of the super heavy neutrino corresponds to the values of the new mass parameters in 
the range St <C Sw *C 1. 

One can prove by means of the direct calculation that 



lim 5w3 n = 0, and lim 8i3 n =0, n = 0, . . . , 2, (4.20) 

5w->0 <5 f — o 



where 



dzz n — . (4.21) 

Using Eqs. 1(4.20(1 and ((4.211) . we find that the function jl in Eq. 1(4.19(1 is equal to 1/2 that corresponds to the magnetic 
moment 

M=^m, (4.22) 



Eq. 14.22(1 presents the magnetic moment of a heavy neutrino with the mass much greater than W boson mass. 



19 



At the end of this section let us compare the calculation of the neutrino magnetic moment in the unitary and i?j 
gauges. The calculations performed within these gauges, as it was mentioned in Ref. |20j, are formally equivalent, 
i.e. the two Feynman amplitudes become equal if we approach the limit a — > oo prior corresponding loop integrals 
are carried out. The diagrams involving unphysical scalar bosons must disappear in the unitary gauge. Therefore, 
in such diagrams the limit a — > oo and the integration over virtual momenta must be commuting procedures. We 
directly verified this statement for particular case of the calculation of the massive neutrino magnetic moment. Indeed, 
on the basis of either exact formulae (|4.3f) - lfi~7|l or the expan sion s giv en by Eqs. H4.9fH4.13f) . we find out that the 



contributions of the diagrams depicted in FigT [l(b)j |T(cj) [Tfe 
limit a 



and 



1(f) which involve the scalar boson, vanish in the 



oo. 



A. Neutrino Magnetic Form Factor at Non-zero Momentum Transfer 



In this subsection we study the neutrino magnetic form factor at non-zero momentum transfer in arbitrary 
gauge as well as for arbitrary charged lepton and neutrino mass parameters (a and b) . The 7 — Z self-energy diagrams 
shown in Fig. 2(a)|2(h) also do not contribute to the magnetic form factor at q 2 7^ 0. Therefore, the total one-loop 
value for the neutrino magnetic form factor is given by 



/m(<7 2 ) 



eG F 



4ir 



v i— 1 



where fj$ (q 2 ) are the contributions to the magnetic moment from the corresponding diagrams shown in Fig. l(a)|l(f) 



For the coefficients fj$ (q 2 ) we have 



[ dz [ dy{2 - 3z + z 2 )^- — 
Jo Jo D i 

dy{az 2 - bz 2 {l - z) - ty(z - y){2 - z)) 



dz 



Di(a) Di_ 

- [ dz [ dy{2 - 3z) [lnDi(a) - hxD x ] , 
2 Jo Jo 



(4.23) 



/m(9 2 ) = ^/ dz I dyz{a + az-b{\-z))-^- 
2 Jo Jo £>iW 



(4.24) 



1 / f Z 1 

f$ ( q 2 ) = - dz dy(2a - 3az + az 2 -bz(l- z)) 



Jo 



D 2 { a y 



(4.25) 



/mV) = o / dz I dyz(l + 2z)—+ 



1 



Jo 



1 



D 2 



dz / dy(b(l ~ zY(z(l - z) - 2y) - ty(z - y)(2y - 3z + z 4 ) - 2ty) 



Jo 



_D 2 (a)+y(l-a) D 2 

If dz f dy(-2 + 9z- Az 2 - 6y) [In (D 2 (a) + y(l - a)) - lnD 2 ] - 
2 Jo Jo 

1 1 2 



dz / dy(bz(l - 3z + z z + z A ) - ty(z - y){2 - z - z 1 )) 



D 2 D 2 {a) D 2 (a)+y(l-a) 



dz dy(8- 13^ + 3^) []nD 2 + \nD 2 (a) -2\n.{D 2 {a) +y(l -a))], (4.26) 
Jo 
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/i? +(6 V)= dz dyy 



i rz 1 



o 



o D 2 (a) + y(l - a) 

1 



i 

2 



dz / dy((a - bz)(l - zf + ty(z - y)(l - z)) 

I! JO 



D 2 (a)+y(l-a) D 2 (a) 



1 

2 Jo 



i 



dz dy(2-3z) [\n(D 2 (a) +y(l- a)) -hft(a)], (4.27) 



u 



where 



Z?i(ct) = a + (a — a)z — bz(l — z) + ty(z — y), D\ = D\(a = 1) = 1 + (a — l)z — bz(l — z) + ty(z — y), 
D 2 (a) = a + (a — a)z — bz(l — z) + ty(z — y), D 2 — D 2 (a = 1) = a + (1 — a)z — bz(l — z) + ty(z — y), 

and t = -q 2 /M^. 

We discuss below the large positive t behavior of the integrals in the expressions of the proper vertices contributions 
t° fid (q 2 )- For example, let us consider the following integral at t — > +oo: 

J(t) = t f dy-^— = f dy- 1 (4.28) 

Jo D 2 (a) J Q {y-y 2 ){yi~y) 

where 

D D a . . 

yi = z + - + ••• , y 2 = — ^ + • • • • (4.29) 

Performing the integrations we readily find that 

J(t)->lnt-lnD. (4.30) 

In Eq. (|4.30l) we dropped the terms like 1/i and (In t)/t which are negligible for large positive t. The remaining 
integrals are evaluated in a similar way. Finally, we find that 



/M(t) = ^/ii ) (t)^0, ift^+oo. 

i=l 



The behavior of the magnetic form factor at large negative q 2 , described above, is consistent with the general Weinberg 
theorem pll |. However, the case of the massive neutrino magnetic form factor has never been discussed previously. 

It should be noted that in derivation of Eqs. 14.28|l - (|4.30|l we assumed that a < oo. Therefore, our result that 
/jvf(i) at t —>■ +oo is valid in any gauge except the unitary one. The value of /m(£ — * +oo) may not be equal to 
zero if we at first set a = oo and then approach the limit t — > +oo. The analysis of the large negative q 2 behavior of 
magnetic form factor within the Weinberg-Salam model in the unitary gauge is given, for instance, in Ref. po| . 

Using the explicit formulae for the massive neutrino magnetic form factor for arbitrary gauge parameter a 
[Eqs. I|4. 2311 - 14. 27fl ] we present in Fig. |3| the behavior of the function fu{t) f° r different gauges and a wide range 
of t: < t < 5 x 10~ 4 . It can be seen that the magnetic form factor becomes gauge independent at t = that 
amounts to the case of on-shell photon. The value fi\i(t — 0) is equal to the neutrino magnetic moment, and Fig. [31 
shows gauge independence of this quantity in agreement with our exact calculations performed above. 



V. CONCLUSION 



We have considered the massive neutrino electric charge and magnetic moment within the context of the standard 
model supplied with SU(2)-singlet right-handed neutrino in general gauge. Using the dimensional-regularization 
scheme, we have calculated the one-loop contributions to the neutrino electromagnetic vertex function taking exactly 
into account the neutrino mass. We have presented the results of our calculations of different contributions to the 
neutrino electric charge and magnetic moment as the closed integral expressions. It has allowed us to determine 
the dependencies of these contributions on the neutrino and corresponding charged lepton masses as well as on the 
gauge-fixing parameters. The integral expressions for the neutrino electric charge and magnetic moment obtained 
in this work contain at most two definite integrals which, in principle, can be performed and expressed in terms 
of elementary functions. However, the results are quite cumbersome and therefore we have presented them as the 
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1.4992 J ■ . . , ^ 

1 2 3 4 5 

**10~ 4 

FIG. 3: The massive neutrino magnetic form factor versus squared transverse momentum in different gauges. The dashed line 
corresponds to a = 100, the solid line to the t'Hooft-Feynman gauge (a = 1), and the dash-dotted line to a = 0.1 

expansion over the neutrino mass parameter b. For several diagrams, which contribute to the neutrino charge and 
magnetic moment and which have been calculated in Ref. 14] with mistakes, we have found the correct results. 

We have found the general expressions for the contributions to neutrino electric form factor. These formulae have 
been derived in general gauge and at arbitrary value of q 2 . We have shown that the electric charge of a massive 
neutrino is a gauge independent and vanishing parameter in the first two orders of the expansion over the neutrino 
mass parameter b. In the particular choice of the 't Hooft-Feynman gauge we have also demonstrated that the neutrino 
charge is zero in all orders of expansion over b, i.e. for arbitrary mass of neutrino. In the previously published works 
devoted to the calculation of the neutrino electric charge the case of massless neutrino was studied within the Georgi- 
Glashow (see Ref. 20]) and Weinberg-Salam (see Refs. 0,0,0,01) models. However, it is clear that the massless 
particle must be electrically neutral. Although there is no doubt that the massive neutrino electric charge must be 
also zero, however it has not been yet shown how it actually happens for corresponding Feynman diagrams. 

There are other reasons to prove by the direct calculations that the value of the massive neutrino electric charge 
is zero. For exam ple, this problem is important in consideration of the neutrino spin oscillations. In the series of 
our works jl^, US EJ we have elaborated the quasi-classical approach for the description of the neutrino spin 
oscillations in arbitrary external electromagnetic field. An essential point in that studies has been the zero charge of 
the massive neutrino. In this paper we have substantiated this assumption. 

From the obtained closed two-integral expression for the massive neutrino electric form factor it is also possible to 
derive the neutrino charge radius. 

The structure of the massive neutrino electromagnetic vertex function have been examined in this work. We have 
directly verified the decomposition of the neutrino vertex function. It has been found out that the value of the 
additional "form factor" fsiq 2 ), which is proportional to 7^75 matrix, at q 2 = is zero in the first two orders of the 
expansion over the neutrino mass parameter b and for arbitrary gauge parameter a. The vanishing of the additional 
"form factor" /5(g ) in the particular gauge (az = 00 and ay? = 1) but for the arbitrary charged lepton a and 
neutrino b mass parameters as well as for arbitrary q 2 has also been demonstrated. Such direct calculations have 
never been carried out before. 

For each of the diagrams contributing to the neutrino magnetic moment, we have obtained the expressions ac- 
counting for the leading (zeroth) and next-to-leading (first) orders in the neutrino mass parameter b with the gauge 
dependence shown explicitly. Each of the contributions is finite and the sum of all contributions turns out to be 
gauge- independent. Our calculations also enable us to get the neutrino magnetic moment in the following ranges of 
the neutrino, charged lepton, and W boson masses: m v <C mi <C My?, me <C m„ <C My?, and mi <C My? <C m v , 
which span almost all the cases presently discussed within different theoretical models. We have also presented the 
general formulae for the massive neutrino magnetic form factor at arbitrary q 2 . 

As for the behavior of the neutrino magnetic form factor at q 2 y^0, we have found that the function /m(<? 2 ) essentially 
depends on the gauge fixing parameter a at q 2 ^0. The magnetic form factor may depend on the gauge parameter at 
q 2 y^0 since it is not a measurable property of a particle and, therefore, may not be invariant under the gauge group 
transformations. The consideration of the gauge parameter dependence of the neutrino magnetic form factor as well 
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as its asymptotic behavior at large negative q 2 in the limit b — > within the Weinberg-Salam model is presented in 
Ref. |20|. The transition magnetic moment of the Dirac neutrinos coupled with the light fermion / (my m m u ) and 
with the light scalar boson <f> (m^ <C m„ — rrif) through the Yukawa interaction Df(f> is discussed in Ref. 26]. The 
transition magnetic moment dependence on q 2 is also considered there. The analysis of the neutrino magnetic moment 
is presented in Ref. [27| for various versions of the left-right symmetric models. The results of our massive neutrino 
magnetic moment calculations can be applied to the treatment of the magnetic moment (including the transitional 
magnetic moment) within the left-right symmetric model. 

Although we have not studied the neutrino anapole form factor in this paper, this particular problem (which we 
discuss in [? ]) is also important since, for instance, the anapole moment (the value of the anapole form factor at 
q 2 = 0) is the only static electromagnetic property of a Majorana neutrino (see, for instance, Refs. |3.l2qV It should 
be noted that even a massless particle can posses the anapole moment, unlike the magnetic moment. Some resent 
papers are worth mentioning in this respect (see plj |29| and references therein). However, the investigation of the 
neutrino anapole moment faces serious difficulties such as its observability and gauge dependence. 
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APPENDIX A: FEYNMAN RULES 



In this Appendix we present the full list of the Feynman rules |3(| necessary for the calculation of the massive 
neutrino electromagnetic vertex. In the gauge the propagators for the vector bosons, W and Z, an unphysical 
charged scalar boson, x, as well as charged ghosts, c and c, are presented in the following form 



- (1 - a) 



k 2 



k 2 



M z z + it 



9fiv - (1 - Q-z) 



aM 2 , 
k^ky 



k 2 - a z M 



Dix) ( k ) = ^T2 — vr— > D^\k) = D^(k) = —- 2 — -g— j- . 

aMyy — k £ — le aikz^ — k z — it 



The fermion propagator has the standard form 

fi + m n 



S(k) 



it 



where n denotes the type of a fermion. 

All vertices can be divided into several classes. We append below the corresponding graphs and Feynman rules for 
each of these classes. 
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(a) g cos e w {{k-p)"tg a P+ 
(p-9)V 7 + (g-fc)V Q } 




5 



(b) e{(fc -p) 7 g Q ' 3 + 

(p-q)V 7 + (<?-fc)V a } 



FIG. 4: |(a)|(b)| vector boson triplex vertices. 




(a) eg cos 0yK X 



FIG. 5: |(a)|vector boson quadruple vertex. 
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(a) (g/V2)ri 



(b) (g/V2ht 




(c) eQ n -y^, n = i, I 




0j 




(d) (g/2cos9 W /)7 a x 
(1/2 -2Q/ sin 2 6> w + 7s/2) 



(e) -(g/2cos6» w )7aX 
(l/2+2Qj sin 2 6W+75/2) 



FIG. 6: |(a)|(e)| one vector boson and two fermions vertices. 



A, 



x- 1 



(a) -(g/2cos6 w )x 
cos 29 w (p - q)a 



(b) -e(p-q) a 



FIG. 7: |(a)|(b)| one vector boson and two scalar bosons vertices. 
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w: 




w- 



(c) -ig sin 2 9 w M z g a f3 



(d) ieM w g a/ 3 



FIG. 8: |(a)|(d)| two vector bosons and one scalar boson vertices. 



X" X 




(rriiP R - miP L ) (miP R - m,iP L ) 



FIG. 9: |(a)|(b)|one scalar boson and two fermions vertices. 



2G 




c ■ c 1 



(b) g cos 9 w p a 



■ • c 




(c) -ep a 

FIG. 10: |(a)|(d)| one vector boson and two charged ghosts vertices. 

X + X~ 




(a) egg a(1 Y. 

COS 29yy I COS dyv 



FIG. 11: |(a)|two vector bosons and two scalar bosons vertex. 



All the momenta of particles associated with vertices are taken to flow in. Qij represent electromagnetic charges 
of the fields tpij in the units of e. ipij representing the three generations of leptons and quarks correspond to usual 
"up" (all types of neutrinos as well as u, c and t quarks; I3 = +1/2) and "down" (all types of leptons as well as d, s 
and b quarks; I3 = —1/2) components of an isodoublet, respectively, I3 is the third component of the isospin. 

The arrow on a line indicates the direction of the flow of a certain quantum number: the charge for W^, x^i the 
fermion number for ip, the ghost number for c, c. The symbol or at the charged ghost lines stands for the sign 
of the charge carried by the arrow. 

APPENDIX B: FEYNMAN INTEGRALS 

In our calculation of Feynman integrals over virtual momenta we use dimensional-regularization scheme with the 
following natural properties of 7-matrix algebra: 

75} = 0, {7 M , 7„} = 2g^, g^g^ = N, 

where N — 4 — 2e is the number of dimensions. 
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The dimensional regularization of the loop integrals in the Euclidian space is performed in the following way: 



{2-kY 



d 4 k 



(27T) 



N 



d N k = 



A 2e 

N 



(27T) 



dn I ^^dk, 



Q(N) 



where fl(N) = 2n N / 2 /T(N/2) is the area of a unit sphere in N dimensions. The dependence of an arbitrary positive 
parameter A, which has the mass dimensionality, is introduced to provide the total dimensionality of an integral. The 
general technique for calculation of various loop integrals in the dimensional regularization scheme can be found, for 
example, in Ref. [3j. It should be, however, rather helpful to include here some of the typical loop integrals which 
one encounters while calculating the electromagnetic vertex function, 
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where 



0.5772157 is the Euler constant. 
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